Accessibility of the Pre-Big Bang Models to LIGO 



Vuk Mandic 

LIGO Laboratory, California Institute of Technology, MS 18-34, Pasadena, CA 91125 



o 
o 



Alessandra Buonanno a ' b 
a Department of Physics, University of Maryland, College Park, MD 20742 
b AstroP articule et Cosmologie (APC) * 11, place Marcelin Berthelot, 75005 Paris, France 

(Dated: February 5, 2008) 

The recent search for a stochastic background of gravitational waves with LIGO interferometers 
has produced a new upper bound on the amplitude of this background in the 100 Hz region. We 
investigate the implications of the current and future LIGO results on pre-Big-Bang models of the 
early Universe, determining the exclusion regions in the parameter space of the minimal pre-Big Bang 
scenario. Although the current LIGO reach is still weaker than the indirect bound from Big-Bang 
nucleosynthesis, future runs by LIGO, in the coming year, and by Advanced LIGO (~2009) should 
further constrain the parameter space, and in some parts surpass the Big-Bang nucleosynthesis 
bound. It will be more difficult to constrain the parameter space in non-minimal pre-Big-Bang 
models, which are characterized by multiple cosmological phases in the yet not well understood 
stringy phase, and where the higher-order curvature and/or quantum- loop corrections in the string 
effective action should be included. 
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I. INTRODUCTION 



The Laser Interferometer Gravitational-wave Observa- 
tory (LIGO) has built three multi- kilometer interferom- 
eters, designed to search for gravitational waves (GWs). 
One of the possible targets of such a search is the stochas- 
tic background of gravitational waves. Many possible 
sources of such a background have been proposed (see, 
e -g-j [J 01 f° r reviews). Some of these sources are 
astrophysical in nature, such as rotating neutron stars, 
supernovae or low-mass X-ray binaries. Others are cos- 
mological, such as the amplification of quantum vacuum 
fluctuations duringinflation 0, Q , phase transitions Q , 
and cosmic strings Q ■ Most of these sources are expected 
to be very weak and below the sensitivity of the LIGO 
interferometers. Furthermore, they are constrained by 
several observations. 

The measurement of the cosmic microwave background 
by the Cosmic Background Explorer (COBE) bounds 
the logarithmic spectrum of gravitational waves |3^| to 
^gw(/)/i?oo < 8 x 10- 14 at - 10- 16 Hz 0, where 
h 100 = iJ /(100km/s/Mpc) w 0.72 is the "reduced" 
Hubble parameter Since in standard (slow-roll) in- 
flationary models, the spectrum produced by the para- 
metric amplification of quantum-vacuum fluctuations |j| 
is expected to be (almost) flat at higher frequencies [l(| , 
a similar bound applies at higher frequencies, as well. In 
some inflationary models in which a cosmological phase 
with equation of state stiffer than radiation comes be- 
fore the radiation era, the spectrum at high frequency 



could increase as function of frequency, thus avoiding the 
COBE bound. For example this happens in quintessen- 
tial inflation ^lj. The GW spectrum could mildly in- 
crease as function of frequency in scenarios in which in- 
flation occurs with an equation of state w < — 1 - 
some examples are given in Ref. [lflj where inflation is 
obtained from a non-canonical Lagrangian. In other sce- 
narios of superstring cosmology, as the cyclic/ekpyrotic 
models [lj], the GW spectrum also increases as function 
of frequency, but its normalization makes it unobservable 
by ground- and space-based detectors. 

The arrival times of the millisecond pulsars can be used 
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to place a bound at ~ 10~ 8 Hz (Tj]: f2Gw(/)^ioo < 
9.3 x 10~ 8 . Doppler tracking of the Cassini spacecraft 
can be used to arrive at yet another bound, in the 
10~ 6 - 10~ 3 Hz band [l6|: Q G w(f)hf 00 < 0.014. The 
Big-Bang Nucleosynthesis (BBN) model and observa- 
tions can be used to constrain the integral of the GW 
spectrum / n GW hj 00 d(\n f) < 6.3 x 10" 6 0, S E3- 
Finally, the ground-based interferometers and resonant 
bars can probe the spectrum of gravitational waves in 
the band 10 Hz - few kHz. The most recent bound from 
LIGO is ^GW^inn < 4.2 x 10~ 4 for a flat spectrum in the 
69-156 Hz band 

In this paper, we focus on the implications of the re- 
cent LIGO result on pre-Big-Bang (PBB) models [l9| . 
and we investigate their accessibility to future LIGO 
searches. The PBB models predict a stochastic GW spec- 
trum whose amplitude can increase as a function of fre- 
quency in some frequency ranges. Hence, they can avoid 
the bounds due to the CMB, pulsar timing, and Doppler 
tracking, and predict relatively large background in the 
frequency band where LIGO is sensitive. In Sec. II we 
briefly review the GW spectrum in the minimal PBB 
models. In Sec. Ill we discuss the latest result from 
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the LIGO search for the stochastic GW background. In 
Sec. IV we study how the new LIGO results, and the 
expected future results, constrain the free parameters of 
the minimal PBB models. In Sec. V we discuss how 
modifications of the minimal PBB model can affect the 
observability of the stochastic GW background. Finally, 
we conclude in Sec. VI. 



tational waves Q| as evaluated in Ref. [2! 
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II. THE GRAVITATIONAL WAVE SPECTRUM 
IN THE MINIMAL P RE-BIG-BANG MODEL 



b(ji) 



In the PBB scenario (see, e.g., [lja|2jJ,|2J|22j), the ini- 
tial state of the Universe is assumed to be the string per- 
turbative vacuum, where the Universe can be described 
by the low-energy string effective action. The kinetic 
energy of the dilaton field drives the Universe through 
an inflationary evolution (henceforth denoted dilaton in- 
flationary phase), which is an accelerated expansion in 
the string frame, or accelerated contraction (gravita- 
tional collapse) in the (usual) Einstein frame. The space- 
time curvature increases in the dilaton inflationary phase, 
eventually reaching the order of the string scale. At this 
point, the low-energy string effective action is no longer 
an accurate description of the Universe, and higher order 
corrections (higher-curvature and/or quantum-loop cor- 
rections) should be included in the string action. These 
corrections are expected to reduce or stop the growth 
of the curvature, removing the would-be Big-Bang sin- 
gularity. The exact evolution of the Universe in this 
high curvature and/or strong-coupling phase (henceforth 
denoted by stringy phase) is currently not known |19|| . 
The end of the stringy phase is what one could refer to 
as the "Big-Bang" — the Universe's transition into the 
radiation phase, which is then followed by the matter- 
dominated and acceleration-dominated phases. 

Although the transition between the inflationary PBB 
phase and the post-Big-Bang phase is not well under- 
stood, some models have been proposed in the litera- 
ture which can partially describe it. In the following, 
we focus on the model derived in Ref. |23j where, in the 
string frame, the dilaton-inflationary phase is followed by 
a phase of constant curvature with the dilaton field grow- 
ing linearly in time. It is then assumed that at the end 
of this stringy phase the dilaton reaches the present vac- 
uum expectation value and stops. This model has been 
denoted in the literature as the "minimal" PBB model. 
Within this model, the stochastic GW background has 
been evaluated 0, HE US • F° r simplicity, in this paper 
we use the result for the logarithmic spectrum of gravi- 
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where Hq 2 { , J M , and T are the Hankel, Bessel and Gamma 
functions, respectively, -ffioo = 100 km/s/Mpc, and Mpi 
is the Planck mass; f s is the GW frequency redshifted un- 
til today of fluctuations exiting the Hubble radius at the 
time of the transition between the dilaton and the stringy 
phase; /i is a dimensionless free parameter that measures 
the growth of the dilaton during the stringy phase, effec- 
tively determining the slope of the spectrum in the high- 
frequency limit (see below). The low- frequency limit of 
Eq. is given by [2^: 
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lE = 0.5772, 
while the high-frequency limit is [2C 
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The parameter /1 appearing in the above equations is 
the GW frequency redshifted until today of fluctuations 
exiting the Hubble radius when the stringy phase ends. 
This is the largest frequency (smallest scale) for which 
fluctuations are amplified — hence, /1 is also the high- 
frequency cut-off of the GW spectrum. 

Thus, the GW spectrum in the minimal PBB scenario 
increases as / 3 for / <C f s , goes as f 3 ~ 2 ^ for f s <C / <C 
/1, and vanishes exponentially for f > fi- An example 
of such a spectrum is shown in Fig. 3 of Ref. [2(| , and 
we reproduce it in Fig. ^ 
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we can define the following cross-correlation estimator: 



Y = 
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FIG. 1: foioo^Gw(/) vs /, as predicted by the PBB model 
with f s = 100 Hz, /i = 4.3 x 10 10 Hz, and fj, = 1.5. 



Let us now focus on the free parameters of the model. 
The parameter /i is, by definition, limited to positive val- 
ues. We will only consider the case fj, < 1.5 — for /i > 1.5 
the decreasing spectrum would easily violate the existing 
experimental bounds |26|. The parameter f s varies over 
the range < / s < f%. Since the spectrum sharply de- 
creases for frequencies below f s , LIGO's reach for models 
where f s is above the LIGO band quickly diminishes. In 
particular, to avoid the / 3 dependence in the LIGO fre- 
quency band, f s < 30 Hz is necessary. Furthermore, 
Eq. @ shows that in the high-frequency limit the spec- 
trum does not depend on f s . Hence, if f s < 30 Hz, it 
does not matter what it is, as far as the accessibility to 
LIGO is concerned. Finally, the parameter j\ can be 
approximated as [2(| |39j : 



/i ^ 4.3 x 10 




(4) 



where H s is the (constant) Hubble parameter during the 
stringy phase, t\ is the time when the string phase ends, 
and A s is the string length. The values H s ~ 0.15A/p; 
and ti ps A s are the most natural ones [2(| 0j but they 
might vary by an order of magnitude. Since flow(f) ~ 
fi [see Eq. (j3Jl], this variation leads to a very large vari- 
ation in the amplitude of the GW spectrum. Hence, al- 
though the theoretically predicted value for /i is more 
robust than those for f s and fj,, we shall explore the pos - 
sibility of varying fi around its most natural value [4JJ . 



III. SEARCHING FOR STOCHASTIC 
GRAVITATIONAL WAVES WITH LIGO 

The method of searching for stochastic gravitational 
waves with interferometers has been studied by many 
authors [H |H . Foilowin g Allen and Romano ET 
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where St is a finite-time approximation to the Delta func- 
tion, §i and §2 are the Fourier transforms of the strain 
time-series of two interferometers, and Q is the optimal 
filter. Assuming that the detector noise is Gaussian, sta- 
tionary, uncorrelated between the two interferometers, 
and uncorrelated with and much larger than the GW 
signal, the variance of the estimator Y is given by: 
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where Pi(f) are the power spectral densities of the two 
interferometers, and T is the measurement time. Finally, 
it can be shown that the optimal filter can be written in 
the form El: 



Q(f) = N 



PPl(f)P2(f) 



(7) 



where 7(/) is the overlap reduction function (arising from 
the different locations and orientations of the two in- 
terferometers), and is the template spectrum to 
be searched. Assuming the template spectrum O t (/) = 
fia(//100 Hz) Q , the normalization constant N can be 
chosen such that < Y >= fi Q T. 

This analysis procedure was implemented in the re- 
cent analysis of the LIGO data, using the 4 km inter- 
ferometers at Hanford, WA and Livingston, LA, for the 
science run S3 [l^. This analysis yielded the 90% upper 
limit of fin, < 8.4 x 10~ 4 for the flat template spectrum 
^t(f) = S^o- Once Y(f) is estimated for the flat spec- 
trum, one can apply simple scaling by the appropriate 
power law to obtain the estimates for different values 
of a (similar procedure can be followed for an arbitrary 
spectral shape). Fig. |2 shows the 90% UL on Q a as a 
function of the spectral slope a for the S3 run, as well 
as the expected reach for LIGO and for Advanced LIGO. 
Here and in the following by expected LIGO (H1L1 and 
H1H2) we mean LIGO design sensitivity and one year of 
observation, and by Advanced LIGO we assume a sensi- 
tivity 10 times better than the LIGO design and one year 
of observation. [LIGO has started the year-long run at 
design sensitivity in November 2005.] 



IV. SCANNING THE PARAMETER SPACE 

We now study the accessibility of the minimal PBB 
model discussed in Sec. II to the most recent and fu- 
ture runs of LIGO, and to Advanced LIGO. Previous 
investigations, which did not use real data, were done in 
Refs. SUD. 



4 



10 ' 



10 



CI 10 



10 



10" 



S3 H1L1 

Expected H1L1 

Expected H1H2 

- - AdvLIGO H1H2 



FIG. 2: The 90% UL on Q a is shown as a function of the 
spectral slope a for the most recent LIGO result. Expected 
sensitivities of LIGO and of Advanced LIGO are also shown. 



As discussed in Section 2, the amplitude of the GW 
spectrum in the PBB models is proportional to / 3 at 
frequencies below f s . Hence, the sensitivity of LIGO to 
PBB models decreases as f s is increased. To avoid the 
/ 3 dependence of the spectrum in the LIGO frequency 
band, we choose f s = 30 Hz. For such choice of f s , the 
LIGO band falls in the relatively flat part of the GW 
spectrum. We vary f\ by a factor of 10 around the most 
natural value estimated in Eq. |@J (i.e., between 4.3 x 
10 9 and 4.3 x 10 11 ) and we vary fi between 1 and 1.5 
(models with /i < 1 are out of reach of LIGO, as shown 
below). For each point in the fi — f\ plane, we evaluate 
fl a = Hgw(/ = 100 Hz) predicted by the model, and 
we check whether it is excluded by the experimental (or 
future expected) results. We also integrate the predicted 
spectrum and check whether it passes the BBN bound 0, 



fi GW (/)/4 d(ln/) < 6.3 x IO- 
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assuming the number of neutrino species N v < 3.9 
[2(||32. We use hwo — 0.72 as the reduced Hubble pa- 
rameter 0. Fig. |3| shows the 90% UL exclusion curves 
obtained in this way. The latest result from LIGO (S3 
run) is just beginning to probe this parameter space. The 
future runs of LIGO (and of Advanced LIGO) are ex- 
pected to probe a more significant part of the parameter 
space, becoming comparable to or even surpassing the 
BBN bound. As expected, LIGO is most sensitive to 
models with fi = 1.5, which corresponds to the flat spec- 
trum at high frequencies. As fi decreases from 1.5, the 
spectral slope increases, and the spectrum in the LIGO 
band drops quickly below LIGO sensitivity. Although 
the BBN bound also weakens for fi < 1.5, the effect is 
not as dramatic because this bound is placed on the inte- 
gral of the spectrum over a large frequency range. Note 
that the LIGO S3 run is sensitive to PBB models with 



fi > 2.7 x 10 11 Hz, relatively large compared to the most 
natural value estimated in Eq. This is true inde- 

pendent of f s : for f s < 30 Hz, Fig. |21 would not change, 
while for f s > 30 Hz all bounds would weaken. Finally, 
the Advanced LIGO is expected to reach models with the 
most natural value of /i = 4.3 x 10 10 Hz. 
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FIG. 3: The 90% UL exclusion curves are shown in the fi — ji 
plane for f a — 30 Hz (the excluded regions are above the cor- 
responding curves). We show the latest result from LIGO, and 
the future expected reach of LIGO and of Advanced LIGO. 
The limit from the BBN is also shown. The horizontal gray 
dashed line denotes the most natural value of /i , given by Eq. 
0. 

It is also possible to use Eq. to turn a bound on 
/i into an exclusion curve in the ii/A s vs H s / (Q.15Mpi) 
plane. In this way, the GW experiments can be used 
to constrain string-related parameters in the framework 
of the PBB model. As an example, we choose fi = 1.5 
and f s = 30 Hz as the optimal case for LIGO, and de- 
termine the 90% UL exclusion curves for different ex- 
periments. These curves are shown in Fig. 0] Again, 
the latest LIGO result is weaker than the BBN bound, 
but the future LIGO and Advanced LIGO searches are 
expected to explore a larger, more physical part of this 
parameter space. 

One can also examine the accessibility of the models 
in the f s — fi plane. For the relatively large value f\ = 
4.3 x 10 11 Hz, which makes the model's stochastic GW 
background accessible to the LIGO S3 run, we performed 
a scan in the / s — fi plane. Fig.[S]shows the results. Note 
that for the flat spectrum (/x = 1.5), the S3 run of LIGO 
is sensitive to models with f s < 120 Hz; future runs of 
LIGO and Advanced LIGO are expected to probe higher 
values of f s as well. Also note that the exclusion curves 
in Fig. [3 are almost vertical (i.e. not very sensitive to 
f s ). This is a consequence of the large value of f\ - for 
smaller values of f\ , the accessibility of models to LIGO 
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FIG. 4: The 90% UL exclusion curves are shown in the ii/A s 
vs H s /(0.15Mpi) plane, for n — 1.5 and f s = 30 Hz (the 
excluded regions are above the corresponding curves). We 
show the latest result from LIGO, and the future expected 
reach of LIGO and of Advanced LIGO. The limit from the 
BBN is also shown. The black circle denotes the most natural 
point, as given in Eq. 
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FIG. 5: The 90% UL exclusion curves are shown in the f a —\i 
plane, for /i = 4.3 x 10 11 Hz (the excluded regions are to 
the right from the corresponding curves). We show the latest 
result from LIGO, and the future expected reach of LIGO and 
of Advanced LIGO. The indirect limit from the BBN excludes 
the whole region shown in this plane. 



would depend more strongly on the value of / s . 

Several papers in the literature 0i HE E3- 
parametrize the GW spectrum in the minimal PBB 
model in terms of z s = fx/ f 8 and g s , defined by g 8 /g\ — 
(fs/fi) , with f3 given by 2^ = 1 2/3 - 3|. The parame- 
ter z s is the total redshift during the stringy phase, thus 
it quantifies its duration, while g\ and g s are the string 
couplings at the end and at the beginning of the stringy 
phase, respectively. Fig. shows the curves from Fig. 
converted into the z s — g s plane, using fi = 4.3 x 10 11 
Hz, and setting g\ to its most natural value given by 
5i/(4tt) = "gut- 



V. GOING BEYOND THE MINIMAL 
P RE-BIG-BANG MODEL 

In this section we investigate how extensions of the 
minimal PBB model or variations of it can impact the 
accessibility of the stochastic GW background to LIGO 
and to Advanced LIGO. 

The GW spectrum in the minimal PBB model was 
originally evaluated |24l 12a . |26j neglecting the higher- 
curvature corrections in the equation of tensorial fluctu- 
ations during the stringy phase. Gasperini [34| evaluated 
the higher-order equation for tensorial fluctuations and 
showed that these corrections modify the amplitude of 
the perturbation only by a factor of order one. Hence, 
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FIG. 6: The 90% UL exclusion curves are shown in the z a —g s 
plane, for /i = 4.3 x 10 Hz. We show the latest result 
from LIGO (thick solid), and the future expected reach of 
LIGO (thin solid for the H1L1 pair, dashed for the H1H2 
pair) and of Advanced LIGO (dash-dotted). The two sets of 
curves correspond to positive (left) and negative (right) signs 
of (2/3 -3). 
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these corrections are not expected to affect our results 
significantly. 

In Refs. [33L [35| the authors have examined the effect 
of radiation production via some reheating process occur- 
ing below the string scale. Such a process may be needed 
to dilute several relic particles produced during (or at 
the very end of) the PBB phase, whose abundance could 
spoil the BBN predictions [36J. Depending on when and 
for how long the entropy production occurs, it can change 
both the shape and the amplitude of the GW spectrum in 
the frequency region around 100 Hz. In general, the am- 
plitude of the spectrum at these frequencies is reduced. 
If we assume that the reheating process occurs at the end 
of the stringy phase (i.e., all of the entropy is produced 
at the end of the stringy phase), then the effect of the 
process is a simple scaling of the original spectrum by the 
factor (1 — 5s) 4 / 3 , where Ss is the fraction of the present 
thermal entropy density that was produced in the pro- 
cess. Fig.[7|shows the exclusion curves in the f s — [i plane 
for Ss = 0.5. By comparing to Fig. O we can see that the 
effect weakens all bounds. 
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FIG. 7: The 90% UL exclusion curves are shown in the f s — jj, 
plane for /i = 4.3 x 10 11 Hz and for 8s — 0.5 (the excluded 
regions are to the right of the corresponding curves). We 
show the latest result from LIGO, and the future expected 
reach of LIGO and of Advanced LIGO. The indirect BBN 
limit excludes all models shown in this plane. 



spectrum of this model can, therefore, be approximated 
by: 



«GW(/) = 



r> DO (£) c 
o 



f<fs 
f>fs 



(9) 



The BBN bound becomes weaker because the integral in 
Eq. (jSJ) is performed over a much smaller frequency range, 
and it can be written as Hdo < 3.8 x 10~ 5 . Fig. [S]shows 
the bound from the latest LIGO result as a function of 
f s . Note that this bound is already better than the BBN 
bound for / > 300 Hz. 




100 



200 



300 
f (Hz) 



400 



500 



FIG. 8: The 90% UL on Q.do is shown as a function of f s 
for the models where stochastic gravitational background is 
not produced during the string phase. The latest LIGO result 
and the BBN bound are shown. 

Finally, as first noticed in Ref. 35], it is well possible 
that many more cosmological phases are present between 
the pre- and the post-Big-Bang eras — some examples 
are given in Refs. [35l l37|. If this is the case, the GW 
spectra during the high-curvature and/or strong coupling 
region will be characterized by several branches with in- 
creasing and decreasing slopes. Due to the dependence of 
the spectra on a larger number of parameters, it would be 
more difficult to constrain these non-minimal scenarios, 
even when LIGO overcomes the BBN bound. 



Another possible, but somewhat arbitrary variation of 
the model, was examined by Allen and Brustein |30| . 
They assumed that stochastic gravitational waves are 
not produced during the stringy phase, but only during 
the dilaton phase. This is achieved by setting /i = / s 
and assuming that Qqw vanishes for f s < f. Such a 
model is not well motivated in the PBB scenario, but it 
is phenomenologically interesting as it represents a class 
of models whose spectrum peaks in the LIGO band. The 



VI. CONCLUSIONS 

Using the most recent LIGO search for the stochastic 
gravitational background we determined the exclu- 
sion re gion s in the parameter space of the PBB minimal 
model 19]. We found that the most recent S3 run can 
access the stochastic GW background only if /i is larger 
than the most natural value 4.3 x f0 lo Hz (i.e. only if 
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/i > 2.7 x 10 11 Hz; see Fig. 01. In this case the S3 run 
of LIGO can exclude the region in the parame- 
ter space with [i 1.5 (i.e., almost a flat spectrum) and 
f s < 120 Hz (see Fig. |SJ|. A one- year run of LIGO at the 
design sensitivity will be able to start excluding regions 
with sligthly increasing GW spectrum (as a function of 
frequency), while Advanced LIGO could exclude spectra 
with slopes of at most « 0.5. Models with larger values 
of f s will become accessible, as well as with lower values 
of /i, including the most natural value fi = 4.3 x 10 10 
Hz. 

As shown in Fig. [31 the BBN bound already excludes 
all models accessible to the LIGO S3 run. However, it 
should be noted that: (i) the LIGO S3 bound is a result 
of a direct measurement of the stochastic background of 
gravitational waves while BBN bound is not, and (ii) 
future searches by LIGO and by Advanced LIGO are 
expected to approach and even surpass (in some parts of 
the parameter space) the BBN bound. 

Analysis of the search in the parameter space more 
commonly used in the literature (see Fig. |SJ) shows that 
LIGO and Advanced LIGO can bound the duration of the 
stringy phase and the string coupling at the beginning 
(end) of the stringy phase (dilaton inflationary phase). 
Similarly, as shown in Fig. 01 by constraining /i these ex- 
periments can constrain other string-related parameters, 
such as H s (the Hubble parameter during the stringy 



phase) and ti/X s (the ratio of the end-time of the stringy 
phase and of the string length) or the value of the string 
coupling at the end of the stringy phase g\ . 

As emphasized above, the stringy phase is not well 
understood, yet. Many variations to the minimal PBB 
model analyzed in this paper are possible and have been 
proposed |U • They can significantly change the 

shape and the amplitude of the spectrum in the frequency 
range around 100 Hz, hence improving or reducing the 
accessibility of the PBB models to LIGO. The presence of 
multi cosmological phases [35l | during the stringy phase 
will make much harder the determination of the exclu- 
sion regions in the PBB parameter space. More robust 
predictions for the stringy phase would be strongly de- 
sirable. 
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